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We developed a detailed microscopic method that describes the shear modulus anomaly of solid 
helium at low temperature. The shear modulus was calculated using the pinning length of dis- 
locations determined in detail for both crossing network nodes and 3 He impurities. The strong 
suppression of the shear modulus is reproduced well as the temperature or stress increases. The 
shear modulus at low temperatures depends strongly on how the state was prepared. All the key 
features in the stress hysteresis can be quantitatively explained in terms of the thermomechanical 
path-dependent pinning length of dislocation networks. 



I. INTRODUCTION 

The observation of a shear modulus anomaly in solid 
4 He below 200 mK has attracted attention because of 
its remarkable similarities in temperature, 3 He concen- 
tration, and drive dependence to the torsional oscillator 
response, which was previously interpreted as the appear- 
ance of non-classical rotational inertia (NCRI)P^ The 
shear modulus anomaly can provide an alternative ex- 
planation within the classical framework for the resonant 
period drop without invoking superfluidity. 

Despite the remarkable similarities)^' the fundamen- 
tal relationship between the shear modulus and the NCRI 
is not clear. For instance, no evidence of a period drop 
was detected, although shear stiffening was also observed, 
in hexagonal closed packed 3 He at low temperatures.^ 
Namely, decoupling can be observed only in a Bose solid, 
whereas stiffening occurs regardless of the quantum prop- 
erties of solid helium. In addition, the response of the 
NCRI and shear modulus anomaly to DC rotation ap- 
peared to be very different. When DC rotation was su- 
perposed on top of AC oscillation, NCRI was dramati- 
cally suppressed, whereas the shear modulus was unaf- 
fected, indicating that the shear stiffening has a different 
microscopic origin from that of NCRlP 

The characteristic features of the shear modulus 
anomaly can be explained by simple classical mecha- 
nisms without including superfluidity. The elastic prop- 
erties of solids were originally explained by the motion 
of dislocations.^ Dislocations, common line defects in a 
solid, are the source of slippage, which softens a crys- 
talline solid. In case of hep helium, the dislocations are 
easily movab le in the basal plane (0001) of the hexago- 
nal structure! 10 ' 11 ! This dislocation motion leads to the 
reduction of the elastic coefficient C44 and consequently 
affects the temperature-dependent shear modulus. The 
stiffening of a solid at low temperatures can be attributed 
to the suppression of the motion of dislocations, which 
can occur when the dislocation lines are pinned by im- 
purities in the solid. These impurities can be detached 
by large strains or thermal fluctuations, which soften the 
solid. 



When solid helium is not rigid but viscoelastic, its re- 
sponse to external AC agitation can differ from that of 
a rigid solid. The response from a viscoelastic solid is 
characterized by the relaxation model with Debye sus- 
ceptibility, in which the shear modulus and dissipation 
are obtained as the re al and imaginary parts of the sus- 
ceptibility, respectively! 12 ' 13 lThe characteristic relaxation 
time grows at low temperature and causes a crossover 
from the relaxed to the unrelaxed phase. The drive de- 
pendence and hysteresis can then be understood in the 
context of a relaxation response. However, the expected 
change due to dissipation in solid helium is too small 
to explain the large shift in the shear modulus at low 
temperatures.^ A more complicated type of suscepti- 
bility, such as glassy susceptibility with suitable extra 
parameters, is necessary to understand the inconsistency 
between the quanti tative magnitude of the shear modulus 
and the dissipation! 14 * 15 ! The shear modulus and the dis- 
sipation with a relatively reasonable ratio were recently 
predicted by applying a distribution of the binding en- 
ergy between 3 He impurities and dislocations.^ 

Iwasa, on the other hand, suggested a simplified model 
to describe the temperature- and 3 He-concentration- 
dependent shear mo dulu s on the basis of the Granato- 
Lucke (GL) modelP^! In this vibrating dislocation 
model, the vibration of dislocations is essentially char- 
acterized by the average dislocation length, which was 
approximated by the average network pinning length at 
high temperatures and by the average impurity pinning 
length at low temperatures. The crossover is induced 
by progressive shortening or lengthening of the average 
pinning length. This mostly gradual change in the total 
fraction of 3 He impurities binding on the dislocations is 
determined by the temperature and 3 He concentration. 
The shorter pinning length at low temperatures restricts 
the motion of dislocations and, accordingly, stiffens the 
solid helium. The dissipation, which is equivalent to the 
damping of the dislocation motion in this model, shows 
good agreement with the observed experimental results. 

Although all the models mentioned above qualitatively 
explain the characteristics of the shear modulus, the un- 
usual drive- and temperature-dependent hysteresis is dif- 
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FIG. 1. Schematic stress-strain relation for a single network line using (a) the original Granato-Lucke model and (b) the 
modified model, (c) Shear modulus for a single network line calculated using the modified model. 



ficult to understand with these simplified models. Here, 
we investigated the microscopic description of the drive- 
and temperature-dependent hysteresis on the basis of dis- 
location pinning by 3 He impurities. Experimental re- 
sults, including the hysteresis, drive dependence, tem- 
perature dependence of the shear modulus, and dissipa- 
tion, are reproduced entirely with a single set of fitting 
parameters. 



II. THE MODEL 

A. Original Granato-Lucke model 9 

Dislocations in a solid constitute a complicated three- 
dimensional network. The motion of dislocations is 
strongly pinned in the nodes (the intersections of the net- 
work). The dislocations can be more weakly pinned by 
impurities, which depends strongly on the temperature 
and stress in the solid. Accordingly, the dynamics of the 
dislocations can alter the elastic properties of a solid sig- 
nificantly. The external shear stress a applied to a solid 
induces a strain that can be represented by the sum of 
the elastic strain and an additional strain due to the vi- 
brating dislocation loops, 



(-total — (el + Erfi 



(1) 



Because the contribution of a single dislocation loop of 
length L to the strain is proportional to L 3 , the total 
additional strain from dislocations is simply expressed as 
the sum of each contribution^ 
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where C 
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v is Poisson's ratio (= 0.3), R is the 



orientation factor, and fj,Q is the elastic shear modulus. 
For the ramdom orientations in polycrystalline solids, the 
anisotropic effect can be averaged into R (~ 0.2)PlThis 
dislocation-induced strain decreases the shear modulus 
compared to that in the perfect crystalline solid with no 
dislocation. Using Eqs. (1) and (2), one may write the 



normalized shear modulus (the shear modulus /i divided 
by the elastic shear modulus /io) in terms of the loop 
length, 
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In the absence of impurity atoms, the characteristic dis- 
location loop length is determined solely by the length 
L n of the dislocation network link between two network 
nodes. When a sufficient number of impurity atoms, such 
as 3 He atoms, are present, they diffuse into the strained 
region near the dislocation lines and pin the dislocation 
lines, functioning as pinning centers. Consequently, the 
average length of the dislocation loops is reduced and 
gradually approximates the impurity pinning length L^ 
determined by the loop length between two impurity pin- 
ning sites. Because the additional strain is directly re- 
lated to the length of the dislocation loops, the pinning 
of dislocations by impurities will suppress the total strain 
significantly. Accordingly, the shear modulus of a solid 
with fully pinned dislocations is nearly the same as that 
of a perfect crystal. 

Given that the external shear stress on a pinning point 
is greater than the binding force between a dislocation 
loop and an impurity, the stress can break the disloca- 
tion away from the impurity pinning point. Because the 
tension on a pinning point induced by shear stress is pro- 
portional to the total length of two neighboring loops, 
Ln + Li2, dislocations shorter than a characteristic criti- 
cal length L c remain pinned. The critical length is deter- 
mined by the maximum binding force / m , and external 
stress a, given by L c (a) = 7r/ m /4&(7, where b is the mag- 
nitude of Burgers vector. The breakaway condition in 
terms of the critical length is, therefore, 
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(4) 



Granato and Lucke pointed out that this stress-induced 
breakaway is a catastrophic process. For instance, once 
the breakaway condition is satisfied, a dislocation will be 
detached from the corresponding pinning site and pro- 
duce a longer loop (Ln + L i2 ). The resultant loop will 
constitute a new longer pair with another adjacent loop 
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(Ln + Li2,L i3 ) that obviously satisfies the breakaway 
condition, Ln + Li2 + L^ > L c . This chain-reaction pro- 
cess continues catastrophically until the entire impurity- 
pinned dislocation is unpinned. Because any pair of ad- 
jacent loops that comply with the breakaway criterion 
can trigger this catastrophic development, the sufficient 
condition for this transition in a dislocation line is given 
by 
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where the subscript max indicates the maximum resul- 
tant length among all the adjoining pairs in the disloca- 
tion line. 

On the other hand, a gradual decrease in stress causes 
progressive pinning of dislocations. The stress-strain re- 
lation for a single network link in the GL model is delin- 
eated in Fig. [ija). 



B. Modified Granato— Lucke model for solid helium 

The GL model is a pioneering model explaining the 
macroscopic physical quantities using the microscopic be- 
havior of dislocations. Nevertheless, minor modification 
is unavoidable when the model is applied to the solid 
helium system. As shown in Fig. [lja), dislocations 
are not pinned until the applied stress decreases back 
to zero in the original GL model. This is because the 
point defects in the solid are not mobile. We assume 
that 3 He atoms, the prevailing impurities in solid 4 He, 
can approach dislocations via quantum diffusion, unlike 
impurities in a normal solidP^ This peculiar property al- 
lows pinning of dislocations even under non-zero stress, 
as shown in Fig. [fjb). The pinning rate can be deter- 
mined by the competition between the binding energy 
of impurities to dislocations and the vibrational energy 
that tears impurities from dislocations. Suppose a 3 He 
atom sticks to a dislocation link that is pinned to two 
neighboring network nodes only. The tension, which is 
greater than the binding force, will remove the 3 He impu- 
rity from the dislocation link immediately after binding. 
The dislocation link can remain pinned only if the ex- 
erted tension is sufficiently small. As described in the 
unpinning process of the original GL model, the tension 
is proportional to the total length of two adjoining im- 
purity pinning loops, which in this case is the length of 
the network link. Hence, the pinning condition can be 
proposed in terms of the critical length as well: 



(6) 



Once a dislocation network link is pinned by a 3 He 
atom, the link will be separated into two shorter loops, 
Ln and Li2, both of which again satisfy the pinning con- 
dition, Ln,Li2 < L c . Therefore, the pinning process in 
solid helium is also catastrophic. For easy comparison, 
the stress-strain relation for a single network line using 
the modified GL model is delineated in Fig. [IJb). 



In Fig. [ljc), the shear modulus [i is obtained as a func- 
tion of the stress by simply using the relation, [i = a/e, 
where a is the external stress, and e is the total induced 
strain. Before the breakaway point in a stress up-scan, 
a single dislocation link shows a relatively high, stress- 
independent shear modulus As the stress reaches 
the breakaway condition, the dislocation link is softened 
by the unpinning of the 3 He impurities, with a lower 
shear modulus ([± u )- The abrupt change of the shear 
modulus in a stress down-scan can be explained in es- 
sentially the same way. As the stress decreases below 
the pinning point, the dislocation is stiffened again ow- 
ing to the reduced loop lengths and recovers the previ- 
ous value of the shear modulus (fi p ). Note that the hys- 
teretic behavior of the shear modulus described in Fig. 
[IJc) is not completely new but rather is similar to the 
hysteretic el astic unit (HEU) in the Preisach-Mayergoyz 
space modelP^E^ In this model, the entire system is as- 
sumed to consist of a large number of HEUs, and its 
macroscopic elastic properties are considered to be the 
integral response of individual HEUspSl 



III. DISTRIBUTION OF DISLOCATION 
PINNING LENGTH 

The hysteresis of the shear modulus shown in Fig. [IJc) 
is limited to a model with a single dislocation link. A 
more realistic description should consider the distribu- 
tions of both impurity pinning lengths Li and network 
pinning lengths L n . The loop length distribution dra- 
matically changes the shape of the hysteresis, as shown 
in Fig. g 

For randomly arranged solvent and impurity atoms 
along a one-dimensional dislocation line, the distribution 
of the length between the impurities (which is equivalent 
to Li) follows an exponential functionPS The distribu- 
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FIG. 2. Distribution effect on the shear modulus as a func- 
tion of stress. 
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tion function of the impurity pinning loops is written 
as N(Li)dLi = ^4- exp (— where A is the to- 

tal length of dislocations per unit volume, i.e., the dis- 
location density, and Li a is the average length of Li. 
When network pinning is superposed on impurity pin- 
ning, the pre-existing impurity pinning loops are further 
separated into shorter loops by network nodes, result- 
ing in a more complicated distribution rather than the 
simple exponential form. Only the average value of the 
pinning length L — £"+1- ^ as hitherto been exploited to 
calculate the temperature- and impurity-concentration- 
dependent shear modulus because of this complexity.^ 

This oversimplified model cannot accurately describe 
the real solid system, as illustrated in Fig. [2j For in- 
stance, the smoothed transition from the pinned to the 
unpinned state with increasing applied stress does not 
appear in this model. Thus, we developed a more realis- 
tic model that allows us to consider various network and 
impurity pinning lengths of dislocation lines. 

First, all dislocation segments between nodes in the 
network can be detached and then connected to recon- 
struct a one-dimensional dislocation line with length A. 
All the network nodes can be superposed on this one- 
dimensional dislocation line. The links between the nodes 
contain the original length distribution of a specific dis- 
location length L n . When the dislocation links are ar- 
ranged in ascending order of length, the links with similar 



length can be categorized in a group. We will consider 
a specific group of dislocation links with a network pin- 
ning length L n initially. Later, the discussion will be 
extended to the entire dislocation network to deal with a 
more realistic model. 

Next, additional pinning of dislocation links by impu- 
rities should be considered. The impurity pinning length 
Li is defined by the dislocation segments pinned by two 
neighboring impurities. For a given specific group of dis- 
location links with a network pinning length L n , we find 
that the effective pinning length cannot exceed L n be- 
cause the impurity pinning length Lj longer than L n will 
be split by network nodes. Specifically, it is divided into 
several segments of dislocation with length L n and seg- 
ments shorter than L n . In addition, Li shorter than L n 
can also be divided by the presence of network nodes 
when a set of impurity distributions is superposed with 
the specific network pinning length L„. (See Appendix.) 
Accordingly, it is necessary to define a new effective dis- 
tribution function for the impurity pinning length in a 
specific dislocation network link, as shown in Fig. [3j 

The new impurity pinning length distribution can be 
obtained by eliminating all the impurity pinning lengths 
split by network nodes from the original distribution and 
adding newly generated dislocation pinning lengths. (See 
Appendix for a detailed derivation.) The new distribu- 
tion function can be written as 



N p (Li, L n )dLi = (1 - j£ + 2J f^)^f- A exp {-^)dL i: for < U < L n 



iV„(L n ) = ^exp(-^) 



(7) 



for Li = L n 



where A n dLi = L n N n (L n )dL n is the total length of the 
dislocations per unit volume (i.e., the dislocation density) 
of which the dislocation links have lengths in the range 
(L n , L n +dL n ); N n (L n )dL n is the distribution function of 
the network links with length L n , which will be discussed 
below, and LiA is the average value of Li . The upper dis- 
tribution function N p {Li, L n ) indicates the effective dis- 
tribution function of impurity-pinned loops with lengths 
in the range (Li, Li + dLi) when the dislocation network 
pinning length is between L n and L n + dL n . The lower 
distribution function, N u (L n ), represents the number of 
impurity-free unpinned links of length L n even at zero 
stress. Note that the distribution function satisfies the 
following normalization condition: 



LiN p (Li, L n )dL t + L n N u (L n ) = A n . (8) 



Now we extend the discussion to the dislocation net- 
work with various dislocation network pinning lengths. 



For a polycrystalline solid, the length of the disloca- 
tion network links are known to follow a log-normal 
distribution,^ 



N n (L n )dL n = Zexp 



(ln L n - ln L r> 

a 2 



dL,- 



(9) 



where s is the width of the distribution, and Z and 
L n are given by - exp(^-) and L n A exp(-^f-), re- 

spectively. L n A is the average value of L n and A = 
J °° A n dL n = f °° L n N n (L n )dL n is the total length of the 
dislocations per unit volume. Ultrasonic measurements 
have revealed that the average network pinning length 
L n A and the width s were several micrometers and in 
the range of (0.5-1.3), respectively.^ Typical values for 
A are 10 9 -10 10 mT 2 for relatively low-quality crystals.^ 
Combining the above distribution functions for Li and L n 
enables us to access the microscopic loop length config- 
urations by performing rather simple integrations. Note 
that the distribution functions satisfy the following nor- 
malization condition: 



5 



A n (L n )dL 

n 



LiN p (Li, L n )dLi + L n N u (L n )]dL n = A. 
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FIG. 3. Effective distribution function due to the superpo- 
sition of the network pinning points on the impurity pinning 
length. 



IV. RESULTS AND DISCUSSION 

A. Effect of temperature and stress on loop length 
configurations 

Conventionally, two types of measurement procedures 
are employed for systematic study of the drive depen- 
dence and temperature dependence of the shear modu- 
lus: (1) a temperature scan and (2) a stress scan. In 
the temperature scan, a cell containing solid helium is 
cooled to the base temperature under a constant stress, 
and then the shear modulus is measured during warmup. 
In the stress scan, the cell is cooled to the desired tem- 
perature under the maximum stress, and then the shear 
modulus is measured during the stress down /up cycle. 
The detailed procedure is described elsewhere P=l 

The microscopic configuration in the pin- 
ning/unpinning process on a network link during 
the temperature scan is depicted in Fig. [lja). At tem- 
peratures higher than approximately 500 mK, 3 He atoms 
evaporate from the dislocations, and the dislocations 
become pinned only by the network nodes or jogsf^ The 
pinning by jogs is not discussed in this manuscript for 
the simplicity. All the pinning loops then have the length 
distribution of network links, N n (L n ), and the resulting 
normalized shear modulus becomes [1 + CRL^\~ l . If 
all the network links have the same constant network 
length L n o, the normalized shear modulus will be 
[1 + Y, CRL lo\~ 1 = I 1 + CRKLl^]- 1 , which leads to 
the previous result, A/x//z = — C-RAL^qP^ Because 
the network nodes function as strong pinning points, 
each L n is considered to be a constant throughout 
the entire experiment unless the dislocation network 
is changed by irreversible plastic deformation. As the 
temperature decreases, because the binding energy 
becomes dominant compared to the 3 He kinetic energy 
of thermal motion, more 3 He atoms tend to pin the 
dislocation lines. As a result, the average impurity 
pinning length LiA decreases depending on the tem- 
perature T according to L^(T) = ^ exp(— ^), where 
a is the interatomic distance in solid helium, cc 3 is 
the 3 He concentration, and Ei, is the binding energy 



between a dislocation link and an impurity atom. The 
decrease in LiA significantly reduces the additional 
strain because the additional strain depends on the 
third power of the dislocation lengths, as mentioned 
above. Consequently, the normalized shear modulus 
increases to [1 + GRL\ (T)] -1 , which is much closer 
to one than the shear modulus at high temperature, 
[1 + ^2,CRL^\~ l . Note that since the relaxation time 
m the 3 He pinning/unpinning process is extremely short 
compared to the measurement frequencies, LiA varies 
continuously as the temperature changes. Therefore, 
no hysteresis appears in the temperature scan, which is 
consistent with experiments.^ 

A similar analogy can be found in the stress-induced 
pinning/unpinning process, as shown in Fig. |4]jb). In the 
high-stress limit, all the links become unpinned, and the 
normalized shear modulus saturates at [1 + Ci?L 3 ] -1 , 
which is independent of temperature. In the low-stress 
limit, all the links become pinned, and the normalized 
shear modulus approaches CRL\ (T)] _1 , which de- 

creases monotonically as the temperature increases. In 
the intermediate stress range, a smooth transition of the 
shear modulus can be found depending on the pinning 
length distributions. Unlike the temperature scan, the 
stress scan always involves the hysteresis in this regime 
owing to the different pinning/unpinning conditions out- 
lined above. Hereafter, the shear modulus for each scan 
will be calculated quantitatively by considering the mi- 
croscopic loop-length configurations. 



B. Pinning process: The stress down-scan 

Assume a dislocation network whose links are of vari- 
ous lengths following a log-normal distribution function, 
N n (L n ). At sufficiently high stress, i.e., for low criti- 
cal length, all the links will vibrate freely without any 
impurity pinning. As the stress decreases, shorter links 
satisfying the pinning condition (L n < L c ) will be pinned 
earlier. When the system reaches the desired stress dur- 
ing the stress down-scan, all the links shorter than the 
critical length will become pinned, but not those longer 
than the critical length. Then the normalized shear mod- 
ulus, Eq. ((3), can be expressed as 
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FIG. 4. Microscopic pinning and unpinning process on a single network link in (a) temperature scan, (b) stress up-scan, and 
(c) stress down-scan. 
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C. Breakaway process: The stress up-scan 

At sufficiently low stress, all the dislocation network 
links will be pinned, satisfying the pinning condition. As 
the stress increases, adjoining loop pairs satisfying the 
breakaway condition ((Ln + Li<i) max > L c ) will trigger 
the breakaway of the entire link including them. The 
probability that a pair (Ln, Li 2 ) stays pinned is the same 
as the probability that the total length of the pair, (Ln + 
Liz), is less than the critical length L c : 



L ° dLn_N v (Ln, L 
Prob p (L n ,L c ) = J ° pK 

r n n IT AT / T 



Jo L py — \ for L n > L c . 

f Q 71 dL a N p (Ln,L n ) x J Q n dL l2 N p (L i2 , L n ) 



(12) 
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The length of any pair on the link cannot exceed L n 
[i.e., (Ln + L i2 ) < L n ]; therefore, when L n < L c , the 
breakaway condition ((Ln + Li 2 ) m ax > L c ) cannot be 
satisfied, and Prob p (L n , L c ) is equal to one. Because any 
adjoining pair on the link could trigger the breakaway, 
the probability that the breakaway is not triggered by 
any adjoining pairs on the link is [Prob p ] k . Here k is 
the average number of adjoining pairs on a link; it is 



given by k = jf^\\ — exp (— 1 , which is the same 

as the average number of impurity points on a link. (See 
Appendix for a detailed derivation.) In other words, the 
[Probp] k fraction of the pinned links stays pinned in the 



Prob p ] k ) fraction 



breakaway process, whereas the (1 — [Prob p ] ) 
becomes unpinned. The loop length distribution for the 
breakaway is then modified into a new form, 



N' p (L h L n )dLi = [Prob p ] k x N p (Li,L n )dL h 



for < Li < L n 



N' u (L n ) = N u (L n ) + (1 - [Prob p ] k ) x _ cxp (_^)] j f or Li = Lr 



(13) 



where ^ [1 — exp(— j^)] is the total number of pinned 
links before the breakaway, which is calculated in Ap- 
pendix. Note that the total length of the dislocations is 
conserved because the reduced length in the upper dis- 



tribution is compensated by the supplemented length in 
the lower distribution. One may find that if L n < L c and 
[Prob p ] k equals one, the loop length distribution reverts 
to the previous distribution form. Then the normalized 
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shear modulus may be expressed as 



Mo 



i+cr( l3 + E ^ 

L n <L c L n >L c 



(14) 



where 
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The normalized shear modulus in the pin- 
ning/breakaway process is given as a function of 
LiAiT) and L c (a). The calculated results are displayed 
in Fig. [5]ja) as a function of T and a instead of LiA 
and L c using the relations, LiA{T) = ^- exp (— ^) 
and -L c (<t) = nf m /4ba, respectively. According to the 
classical elasticity theory proposed by Cottrell, the 
maximum attractive force between an impurity and a 
dislocation is f m = — 3v3£k/8ro, where r is the closest 
distance between the impurity and the dislocation 
center, which is given by b/y/2, and Et is the binding 
energy.^ For the best fit to the experimental data, the 
parameters in the distributions were obtained as follows: 
the binding energy Ej, = 0.3 K, the average network link 
length L n A — 5 /im, the width of the network link length 
distribution s = 1, and RA = 2.01 x 10 10 mT 2 . The 
experimental results are shown in Fig. [5jh) for compari- 
son. The arrows indicate the scan direction during each 
measurement. The stress range and temperature range 
of the crossover from the pinned state with high shear 
modulus to the unpinned state with low shear modulus 
show surprisingly good agreement with the calculated 
results. 



D. Constant stress process: The temperature scan 

The conventional temperature scan resembles the pin- 
ning process in the mechanical breakaway model. Note 
that, in the temperature scan with a fixed stress [Fig. 
|4|a)], the links shorter than the critical length are al- 
ways pinned throughout the experiment, whereas the 
temperature-dependent Li is the only relevant param- 
eter that affects the shear modulus. That is, the disloca- 
tion network link always satisfies the pinning condition 
(L n < L c ). On the basis of this fact, the shear modu- 
lus in the temperature scan procedure is constructed in 
Fig. |6][a) using Eq. (11). The experimental results are 
displayed in Fig. [6jb) for comparison. 

Using the mechanical breakaway model, the tempera- 
ture scan under a constant stress and the stress down- 
scan at a constant temperature always lead to the identi- 
cal configuration and consequently the same shear mod- 
ulus for a certain temperature and stress. The calcu- 
lated results also seem to agree well with the expecta- 
tion. However, we observed a discrepancy between the 
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FIG. 5. (a) Shear modulus as a function of stress calculated 
using Eqs. (11) and (14). (b) Experimental shear modulus 
data measured during the stress down/up-scan! 2 ^ 



shear modulus measured via the conventional tempera- 
ture scan with constant stress and the shear modulus 
measured via the stress down-scan below a certain char- 
acteristic temperature. The deviation provides evidence 
that the low-temperature shear modulus of solid helium is 
hysteretic. The hysteresis between these two procedures 
cannot be understood clearly using only the mechanical 
breakaway model. This issue will be discussed further in 
another article.^ 

The vibrational motion of dislocation is hindered by 
various damping mechanisms, resulting in energy dissi- 
pation. The dissipated energy ratio can be expressed 
as the inverse of the quality factor, which is given by 
Q^ 1 = g^r^r, where AW is the energy lost, and W 
is the total energy of the vibrating loop. According to 
Granato and LuckeP the energy dissipation in terms of 



Q is written as 



Q- 1 = rJ LQ~ 1 {L)N(L)dL. 



(15) 



Here, Q 1 (L) is the energy dissipation due to the dislo- 
cation loop of length L, given by 



Q~\L) 



8^ 



cud 



(16) 



7T 3 /3 (Wq — 0J 2 ) 2 + U! 2 d 2 ' 

where and p are the elastic shear modulus and the den 



sity of solid helium, respectively; ujq is the resonance fre- 
quency of a vibrating dislocation loop (= 



(l-y)p' 

is the frequency of the external stress, and d is the damp- 
ing constant. The primary damping source in solid he- 
lium at sufficiently low temperature is the dragging force 
induced on the dislocation lines by 3 He impurity atomsP^ 
Therefore, the damping constant is proportional to the 
3 He density near the dislocation line, which is equivalent 
to the inverse of the pinning length, 



d{T) = d L. 



(17) 



where do is a fitting parameter, which has been deter- 
mined to be 2.43 x 10 5 . Note that this impurity damping 
occurs only if the dislocation loop interacts with impu- 
rity atoms near the pinning points. When a network link 
is not pinned and vibrates freely, in contrast, the accom- 
panying energy loss would be negligible. On the basis of 
this argument, one may calculate the energy dissipation 
by integrating over all the lengths of pinned loops, i.e., 
the lengths satisfying the condition (Lj < L n < L c ). 



Q- 1 = R 



L 1 Q~ 1 (L i )N p (L i , L n )dLidL n . (18) 



The calculated dissipation as a function of temperature 
is shown in Fig. [7][a). As a higher external stress is 
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FIG. 7. (a) Dissipation as a function of temperature cal- 
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Shear modulus as a function of temperature 
(b) Experimental shear modulus 



(HI 



culated using Eq. (18 1, (b) Experimentally measured phase 
shift during the temperature scanP^ 



applied, the fraction of pinned links, which are the source 
of the dissipation, is reduced, resulting in less dissipation. 
When an external stress is exerted on the solid helium, 
the resulting strain is not in phase owing to the damping 
but reveals some phase lag. The phase angle difference 
(A</>) does not directly indicate the absolute values of the 
energy dissipation but is closely related to the dissipation 
for small phase differences, Q _1 ~ Ac/). The measured 
phase angle difference between the applied stress and the 
measured strain is plotted in Fig. [7|h) for comparison. 



V. CONCLUSION 

We proposed a microscopic description to explain the 
stress- and temperature-dependent shear modulus and 
dissipation in solid 4 He using the mechanical breakaway 
model. A new effective distribution function for the pin- 
ning length involving both network node pinning and im- 
purity pinning was constructed by superposition of net- 
work nodes on the impurity pinning loops. The shear 
modulus is determined by the dynamics of the pinning 
length distribution, which depends on the stress and tem- 
perature trajectory. The hysteresis of the shear modulus 



9 



measured in the stress up/down-scan can be attributed 
to the path-dependent pinning length distribution. Shear 
modulus and dissipation data were successfully repro- 
duced within the framework of the new description and 
exhibited quantitative consistency with experimental re- 
sults. In addition, the effective distribution function of 
the pinning length can provide universal applicability to 
various dislocation network systems. 
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Appendix A: The Effective Distribution 

Assume a subset of dislocation network where all the 
network links have same length, L n . Before considering 
the additional separation of dislocation loops by the net- 
work node, the distribution function for impurity pinning 
length (Li) derived by Koehler can be written as 



N(Li,L n )dLi 



A„ . 
— exp(- 



L 



iA 



J iA 



)dLi, (Al) 



where A n dL n = L n N n (L n )dL n is the total length of the 
dislocations per unit volume and the subscript n indi- 
cates that all the links have the length L n , and L.- l a is 
the average length of Li. The addition of the network 
node splits the impurity pinning loops into smaller loops, 
which changes the shape and range of the distribution 
function of Li significantly. The effective loop length dis- 
tribution modified by the superposition of the network 
node can be discussed for two different cases. 



1. Case 1: For Li < L n 

For Li < L n , several possible configurations of the 
superposition are depicted in Figj8](a). All the config- 
urations in the cycle (A-E) are equally probable since 



(a) Case 1: For i, < L„ 
L„ 



(b) Case 2: Fori, > i„ 
L„ i 



there is no preference for the relative position of the net- 
work node. The addition of a network node between two 
neighboring impurity pinning sites can divide the exist- 
ing dislocation loops into two new dislocation loops as 
shown in the cycle (C-E). The dislocation loops remain 
the same otherwise. Accordingly, the probability to have 
the network node in the dislocation loops is given as the 
ratio of the impurity pinning length, Li, with respect to 
the network link length, L n . Note that the probability 
can be found by considering the possibility for the po- 
sition where the node is on the impurity pinning loop 
(C-E) divided by the possibility for any position in the 
whole cycle (A-E). The new distribution function can 
be obtained by subtracting the eliminated dislocation 
loops from the original distribution and including the 
newly created dislocation loops by the additional sep- 
aration. First, the distribution function change due to 
the elimination, AN- (Li, L n )dLi, is then the product of 
the separation probability and the distribution function 
for the existing impurity pinning loops, which is simply 
Ni(L h L n )dLi: 



AN_(Li, L n )dLi = 



Li 
L n 
U 

L n. 



x N(L l ,L n )dL l 



A„ 



Li 



, exp(-— i-)dL<{A2) 

iA L iA 



Second, the distribution function change by the cre- 
ation of the split dislocation loops can be denoted as 
AN + (Li, L n )Li. Assume an impurity pinning loop of 
length L'i is split into two new impurity pinning loops, 
one of which has the length of Li (apparently shorter than 
L'j). The distribution function for the loops of length L\ 
can be obtained with the same way described above (i.e. 
the product of the separation probability and the distri- 
bution function for the existing impurity pinning loops): 

jf- x N(L i} LnjdLi- Since there is no preference for the 
length of resulting two new loops, the new loops can have 
any length smaller than L[ . The probability that the split 
loop length is found in the range of (Li, Li+dLi) is simply 
given by dLijL\. Hence the distribution function change 
due to the creation can be written by 

rL n jr. tJ 

AN +A (L h L n )dLi = I 2 x — 1 x x N(L' t ,L n )dL[ 



2A„ 



L n L 



iA 



L'i L r 
Li 



L r , 



exp(- — -) - exp(—~) 

J-iiA J^iA 



dL u 
(A3) 



FIG. 8. Schematic configurations of two cases: (a) Li < L n 
and (b) Li > L n . The uppermost line represents the network 
line with network nodes, which is displayed as vertical line 
segments on the line. The horizontal line segments below 
represent the impurity pinning loops in various positions. 



where the factor 2 indicates that always two loops are 
created by the superposition. Note that the range of L • 
is chosen to be (Li, L n ) since only impurity pinning loops 
longer than L t can contribute to the creation of new loops 
of length Li . 
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2. Case 2: For Li > L n 

For the impurity pinning loops longer than L n , the sit- 
uation becomes slightly different. All the impurity pin- 
ning loops should be split by the superposed nodes at 
least once and the resulting loops have length less than 
or equal to L n . Assume an impurity pinning loop of 
length L[(> L n ). The length can be written in the 
following form: 



L\ = mL n + 1, < I < L n , 



m 



1,2,3, 



(A4) 



where to is the number of the resulting loops of length 
L n and I is the loop length of the remainder. 

First, note that two new loops shorter than L n are 
created at both ends of the impurity pinned dislocation 
loop of length L- (see Fig. [f|b)). The length of the 
two split dislocation loops at the both ends varies from 
zero to L n in the (A-E) cycle of which value is equally 
probable between zero and L n . The probability that the 
new loop has a length in the range of (Li, Li + dLi) is 
given as dLi/L n . The distribution function change due 
to the creation of loops at both ends is then 

f°° dL 

AJV + , 2 (Li, L n )dL, = / 2 x -A x N(Li, L n )dL'i 

= J^exp(-^)dL i , (A5) 

where the factor 2 indicates that two loops are created 
from the loop of length L^. Note that the range of L'i is 



chosen to be (L ra , oo) since only the loops longer than L n 
are taken into account here. 

Adding all the distribution function changes to the 
original distribution function, the new distribution func- 
tion is obtained as 

N p (Li, L n )dL l = N(Li,L n )dLi + AN-(Li,L n )dLi 

+AN +tl (Li, L n )dL t + AN + a{L u L n )dL t 

= (1- ^+ 2 -^)^eM-^)dHA6) 

J-'n J^n ^iA ^iA 

Noteworthy is that the new distribution function does 
not have the dislocation loops longer than L n , since the 
pinning length cannot be greater than the network pin- 
ning length. 

However, the dislocation links pinned only with the 
network nodes are not included in the new distribution 
function in the Eq. (A6). The additional creation of 
the loops with L n for an impurity pinning loop of length 
L'i = mL n + I results in either (m — 1) loops of length 
L n (A-C in Fig. gb)) or m loops of length L n (C-E). 
Since the position of nodes can be in any configurations 
in the cycle (A-E), the probability that L^ is divided into 
m loops of length L n is simply the number of possible 
configurations for (C-E) divided by that for (A-E), which 

is given by = Li ™ Ln . The probability for (to — 1) 

loops with length L n is then (1 — 37-) = ( m+1 | L " — ^i. 
Now, the number of resulting loops of length L n created 
by the superposition of nodes on the loop of a length in 
the range of (L^, L^ + dL'j) is calculated by the product of 
the expectation number of new loops and the distribution 
function of the existing loops. 



NL n (L'i, L n )dL'i 



L^ - mL n (to + l)L n - L\ 

TO X h (TO — 1) X — 



= (^--l)^L e xp(-A )dLij 

J^n J^iA ^iA 



N(L'i, L n )dL'i 



(A7) 



Integration of the Eq. (A7) over the range of L'i(L n < 
L'i < 00) will lead to the total number of resulting loops 
of length L n . 



r 



poo A T 

Nu(L n ) = / N Ln (L' z ,L n )dLr = ^exp(-f^) (A8) 

J L n *-"n Li A 



In conclusion, the new effective distribution function is 
written as below. 



N p (Li,L n ) = (1 - ~ + 2 ^)-^cxp(--^)dL l ,for < Li < L n 



N u (L n ) = I ±eM-^ L ) (A9) 

J^n ^iA 



The validity of the effective distribution function can be readily tested with the normalization condition. 
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LiN p (L t , L n )dLi + L n N u (L n ) = A ri 



(A10) 



The first term on the left hand side of the Eq. (A10) is 



the total length of new loops whose length is less than 
L n and the second term is the total length of new loops 
of length L n . 



Mz7i + ^-r^p(-fe)] 

41 - exp(- 



L iA [l-eM-T^)} 



+1 



(B4) 

Since the number of the impurity pinning points in a 
network link is one less than the number of the impu- 
rity pinning loops, the number of the pinning points k is 
therefore 



Appendix B: Characteristics of the Dislocation 
Network 



The new effective distribution function (AlOl is com- 



posed of an impurity-pinned part (Li < L n ) and natu- 
rally unpinned part (L, = L n ) that is only pinned with 
the network nodes. In the following paragraph, the char- 
acteristics of each part will be discussed to help under- 
standing the article. 



L iA [l - exp(-^ 



(B5) 



2. Unpinned part (Li — L n ) 

The total number of the dislocation links of length L n 
with the effective distribution function is 



1. Impurity-pinned part (Li < L n ) 

The total number of the loops shorter than L n with the 
effective distribution function can be obtained by inte- 
grating the distribution function of the impurity-pinned 
part. 



N u (L r , 



A n , L n 

— exp(- — ). 



(B6) 



The total length of the links of length L n with the 
effective distribution function is then obtained simply by 
multiplying L n . 



N p (L i: L n )dL l = A r , 



•1,1 1 / L n 

+ - J- ex P(-^)l 

^iA -^n J-*n ^iA 

(Bl) 

The total length of the loops shorter than L n with the 
effective distribution function can be given by integrating 
the product of the length and the distribution function. 



L n N u (L n ) = A„exp(--^-) 

LiA 



(B7) 



Now one can notice that the total number of the all 
kinds of pinning loops is obtained by adding (Bl) and 
(B6). 



LiN p (Li,L n )dLi = A„[l - exp( 



L 



,A 



-)\ (B2) 



The total number of the impurity-pinned network links 
is then given by the total length of the loops shorter than 
L n divided by L n . 



A n ( 



(B8) 



Moreover, the average length of the pinning loops is 
then 



A n . . L n > 

— [l-exp(— — , 



(B3) 



The average number of the impurity pinning loops in 
a pinned network link may be calculated by dividing the 
total number of the impurity pinning loops (Bl) by the 
number of the pinned links (B3). 



LiALr 



(B9) 



Note that the averaged results, Eq. (B8| and Eq 



(B9), are the same as those obtained in the dislocation- 
vibration modelP^ 
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